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Abstract 



We theoretically study a non-magnetic impurity effect on the vortex bound states 
of a multi-quantum vortex. The zcro-cncrgy peak of the local density of states is 
investigated for vortex cores with the winding numbers 2 and 4 within the framework 
of the quasiclassical theory of superconductivity. We find that the zero-energy peaks, 
which appear away from the vortex center in the clean limit, move towards the 
vortex center with increasing the impurity scattering rate, resolving a contradiction 
between an experimental result and previous theoretical predictions. 
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1 Introduction 

Much attention has been focused on the vortex core structure in superconduc- 
tors. A vortex with a multiple winding number is called the multi-quantum 
vortex or the giant vortex, which accommodates multiple flux quanta. While 
only vortices with the winding number L = 1 are formed in usual situations, 
a multi-quantum vortex with L > 1 may be realized in a small-size super- 
conductor because of a geometrical restriction of the split into individual vor- 
tices with L = 1 [1]. Recently, an observation of multi-quantum vortices by 
scanning tunneling microscopy (STM) was reported, where the experiment 
was performed on nanometer-sized Pb islands [2]. There is, however, a con- 
tradiction between the experimental observation [2] and previous theoretical 
predictions [3,4,5,6,7,8,9]. In the STM experiment, the zero-bias peak of the 
tunneling conductance was observed at the vortex center for both odd and even 
winding- number vortices [2]. In contrast, previous theories predicted that the 
zero-energy peak of the local density of states exists at the vortex center only 
for odd winding- number vortices and not for even ones [3,4,5,6,7,8,9]. There- 
fore, it is necessary to resolve this contradiction for even winding-number 
cases. 

In this paper, we claim that an impurity scattering effect can reconcile the 
above contradiction. The zero-energy local density of states inside multi- 
quantum vortex cores with L = 2 and 4 is investigated under the influence 
of impurities. The zero-energy peaks appear away from the vortex center in 
the clean limit, which is consistent with previous theories [3,4,5,6,7,8,9]. We 
find that these peaks move towards the vortex center with increasing the im- 
purity scattering rate, resulting in a peak situated at the vortex center. The 
experimental result [2] , which seemingly contradicts previous theories, can be 
explained by this impurity scattering effect. 
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2 Formulation 



We consider a single vortex in a single-band s-wave superconductor. The sys- 
tem is assumed to be an isotropic two-dimensional layer perpendicular to the 

vorticity along the z axis. In a circular coordinate system within the layer, the 
real-space position is r = (r cos0, r sin0), and the unit vector k represents 
the sense of the wave number on a Fermi surface assumed to be circular. The 
Fermi velocity is = vpk. 

The pair potential A(r) around a multi-quantum vortex with winding number 
L is represented as 



A(r) = A(r)e^^^, (1) 

where A(r) is the pair potential amplitude and the vortex center is situated 
at r = 0. The vortex core structure is self-consistently calculated by means of 
the quasiclassical theory of superconductivity [10]. The Eilenberger equation 
is numerically solved to obtain the quasiclassical Green's function giibJ^ i", k). 
The effect of impurities distributed randomly in the system is taken into ac- 
count through the impurity self energy S(ia;n, r, k). These quantities and the 
Eilenberger equation to be solved are [10,11,12] 



ivp ■ Vg + [iujnh -A,g]= 0. (3) 

The equation is supplemented by the normalization condition g^ — —tt^To- 
Here, 7^3 is the z component of the Pauli matrix, tq is the unit matrix, and the 
brackets denote the commutator [A, B] — AB — BA. The Eilenberger equation 

contains the renormalized Matsubara frequency a)„ and pair potential A [11], 
•~ • V X ( ^ A + Ei2\ 



We have neglected the vector potential assuming the case of a large Ginzburg- 
Landau parameter. Throughout the paper, we use units in which h — kB — 1- 

Gonsidering an s-wave non-magnetic impurity scattering and the t-matrix, the 
impurity self energy E is given by [12] 
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E(ia;„,r) 



I^^ {-^{g) if) 



) 



(5) 



l-{sm'do){l-C)\-{P) i{g) 



where C = {g^ + (/)(/^) with (■ ■ ■) being the average over the Fermi surface 
with respect to k. The impurity scattering rate in the normal state is F^, 
which is related to the mean free path I — vy/2V^. The scattering phase shift 
is 5q. In the present study, we set 5q — Q keeping finite, which corresponds 
to the Born limit. 

The self-consistency equation for A, called the gap equation, is given as 



where is the cutoff energy and the coupling constant A is given by 



Here, T is the temperature and is the superconducting critical temperature. 
We set Uc — lOAo with Aq being the BCS pair-potential amplitude at zero 
temperature. 

The Eilenberger equation, the impurity self energy, and the gap equation are 
numerically solved self-consistently at T = O.lTc to determine the spatial 
profile of the pair potential amplitude A(r) following a procedure described 
in Rcf. [12]. Then, using the pair potential determined sclf-consistently, the 
Eilenberger equation and the impurity self energy are self-consistently solved 
for a real energy E with replacing iujn E + ir] [13,14]. Here, r/ is a small 
positive energy, which is technically necessary for obtaining numerically a 
retarded Green's function. We set rj — O.OSAq. The local density of states 
A'"(r, E) is calculated as 



where N-p is the density of states at the Fermi level in the normal state. 

In the next section, we will show results for the pair potential and the zero- 
energy local density of states around the multi-quantum vortices with L = 2 
and 4. We define the coherence length = vp/Aq and use it as the unit of 
the length. 




(6) 




(7) 



N{r, E) - NpRe{g{iujn ^E + ir}, r, k)). 



(8) 
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3 Result and Discussion 



In Fig. 1, we show the pair potential amphtudes A(r) around the multi- 
quantum vortices with the winding numbers L — 2 and 4. In the clean limit 
(Fn = 0), the spatial profiles are consistent with results shown in Ref. [6]. A 
healing length, which is a length necessary for the pair potential to recover 
to a balk value, decreases with increasing T^. In other words, the vortex core 
radius shrinks with increasing Fn. 

In Fig. 2, shown are numerical results for the zero-energy local density of states 
A^(r, E = around the multi-quantum vortices with L = 2 and 4. They are 
plotted as a function of the radial distance r from the vortex center for several 
values of the impurity scattering rate Fn- In the clean limit (Fn = 0), the zero- 
energy peaks appear at finite distances of the order of the coherence length, 
which is consistent with previous theories [3,4,5,6,7,8,9]. Wc find that the 
zero-energy peaks gradually move towards the vortex center with increasing 
Fji. The two peaks reduce to a single one in the case of L = 4. In a dirty case 
(Fn = 6Ao), the zero-energy peak is located almost at the vortex center r — 
and the peak height is equal to the normal-state density of states at the Fermi 
level A^F- 

Wc can make a physical interpretation of the result in the dirty state as fol- 
lows. When the quasiparticles inside a vortex core are frequently scattered by 
impurities, they lose the information on the phase of the pair potential that 
they feel. As a consequence, the difference of the winding number L related 
to the pair potential phase becomes ineffective for the electronic structure of 
those quasiparticles and the so-called normal core, as in L = 1, is realized in 
the dirty state irrespective of value of L. (Note that the pair potential itself 
keeps its winding number even in the dirty state, namely the winding number 
is kept being L 7^ and the vortex never disappears.) At the normal core, the 
value of the local density of states is the normal-state density of states A^p- 
On the other hand, far away from a vortex, namely in the bulk, the super- 
conducting gap is fully opened in an s-wave superconductor even in the dirty 
state and the zero-energy density of states is zero. As a result, the zero-energy 
local density of states has the value of Np at the vortex core and decays to 
zero far away from the core irrespective of winding number L. Such a spatial 
variation of the zero-energy local density of states exhibits an appearance of 
a peak structure as seen in Fig. 1 (a dirty case of Fn = 6A0). 

In the STM experiment, the superconducting Pb islands on a Si substrate are 
in the dirty state [2,15]. Therefore, the absence rule of the zero-energy peak 
at the vortex center for multi-quantum vortex with even winding number 
predicted theoretically in the clean limit [3,4,5,6,7,8,9] is not applicable in 
that experiment [2] . The presence of the zero-energy peak at the vortex center 
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for L — 2 observed experimentally [2] can naturally be understood in terms 
of the impurity scattering effect discussed above. 



4 Conclusion 

We investigated the zero-energy density of states around the multi-quantum 
vortices with the winding numbers L = 2 and 4 in the s-wave superconductor 
by means of self-consistent numerical calculations based on the quasiclassical 
theory of superconductivity. We found that the zero-energy peaks, which ap- 
pear away from the vortex center in the clean limit, gradually move towards 
the vortex center with increasing the impurity scattering rate. This result of 
the impurity scattering effect can explain the presence of the zero-energy peak 
at the vortex center for L = 2 observed experimentally [2] . Effects of a finite 
system size [16] and of inclusion of the vector potential are left for future 
studies. 
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Fig. 1. Spatial profiles of the pair potential amplitude A(r) around the multi-quan- 
tum vortices with the winding numbers (a) L = 2 and (b) L = 4 for several values 
of the impurity scattering rate Fn- The horizontal axis r is the distance from the 
vortex center. 



Fig. 2. Spatial profiles of the zero-energy local density of states N{r, E = 0) around 
the multi-quantum vortices with the winding numbers (a) L = 2 and (b) L = 4 for 
several values of the impurity scattering rate Fn. The horizontal axis r is the radial 
distance from the vortex center. 
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Figure 1(a) 



(a) Pair Potential for L=2 
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Figure 1(b) 
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Figure 2(a) 
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(a) Zero-Energy LDOS for L=2 
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Figure 2(b) 
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(b) Zero-Energy LDOS for L=4 
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